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Abstract 

We present a unified derivation of covariant time derivatives, wliicli transform 
as tensors under a time-dependent coordinate change. Such derivatives are essen- 
tial for formulating physical laws in a frame-independent manner. Three specific 
derivatives are described: convective, corotational, and directional. The covariance 
is made explicit by working in arbitrary time-dependent coordinates, instead of re- 
stricting to Eulerian (fixed) or Lagrangian (material) coordinates. The commutator 
of covariant time and space derivatives is interpreted in terms of a time- curvature 
that shares many properties of the Riemann curvature tensor, and reflects nontrivial 
time-dependence of the metric. 



PACS numbers: 83.10.Bb, 05.45.-a, 47.50. +d 
1 Introduction 



In physics, clioosing an appropriate coordinate system can make the difference 
between a tractable problem and one that defies analytical study. In fluid 
dynamics, two main types of coordinates are used, each representing a natural 
setting in which to study fluid motion: the Eulerian coordinates, also known 
as the laboratory frame, are time-independent and fixed in space; in contrast, 
the Lagrangian (or material) coordinates are constructed to move with fluid 
elements. In between these extremes, other types of coordinates are used, such 
as rotating coordinates in geophysical fluid dynamics. Such coordinates usually 
have a nontrivial spatial and temporal dependence. 

The situation becomes more complicated when dealing with moving surfaces: 
here the metric itself has intrinsic time dependence. This time dependence 
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incorporates the strain imposed on a 2D surface flow as the surface deforms. 
To properly formulate fluid equations on thin films and other surfaces, one 
needs a covariant description, that is, a description of the building blocks of 
equations of motion — spatial and temporal derivatives — that obey tensorial 
transformation laws. 

There arc other reasons than inherent deformation of the space to introduce a 
time-dependent, nontrivial metric. For instance, the ad vection- diffusion equa- 
tion can have an anisotropic, time-dependent diffusion tensor, perhaps arising 
from some inhomogeneous turbulent process. In that case, it is advantageous 
to use the diffusion tensor as a metric, for then the characteristic directions of 
stretching, given by the eigenvectors of the metric tensor in Lagrangian coordi- 
nates, correspond to directions of suppressed or enhanced diffusion associated 
with positive or negative Lyapunov exponents, respectively [1,2]. 

Local physical quantities can be viewed as tensors (scalars, vectors, or higher- 
order tensors) evaluated along fluid trajectories. For instance, we may be in- 
terested in how the temperature (scalar) of a fluid element varies along a 
trajectory, or how the magnetic field (vector) associated with a fluid element 
evolves. Characterising the evolution of these tensors in complicated coordi- 
nates is again best done using some form of covariant time derivative, also 
called an objective time derivative. 

The covariant spatial derivative is a familiar tool of differential geometry [3-5] . 
The emphasis is usually on covariance under coordinate transformations of the 
full space-time. In fluid dynamics and general dynamical systems, however, the 
time coordinate is not included in the metric (though the metric components 
may depend on time), and the required covariance is less restrictive: we seek 
covariance under time-dependent transformations of the coordinates, but the 
new time is the same as the old and does not depend on the coordinates. 
Time derivatives lead to non-tensorial terms because of time-dependent basis 
vectors — the same reason that ordinary derivatives are not covariant. 

There are many ways of choosing a covariant time derivative. The most fa- 
miliar is the convective derivative introduced by Oldroyd [6,7] in formulating 
rhcological equations of state. This derivative was then used by Scriven [8] to 
develop a theory of fluid motion on an interface. The convective derivative of 
a tensor is essentially its Lie derivative along the velocity vector. In spite of 
its economical elegance, the convective derivative has drawbacks. Firstly, un- 
like the usual covariant spatial derivative, it is not compatible with the metric 
tensor. A compatible operator vanishes when acting on the metric. Because 
the covariant derivative also has the Leibniz property, compatibility allows 
the raising and lowering of indices "through" the operator. This is convenient 
for some applications [1], and implies that the equation of motion for a con- 
travariant tensor has the same form as the covariant one. A second drawback 



2 



of the convcctive derivative is that it involves gradients of the velocity, and so 
is not directional. The commutator of the convcctive derivative and the spatial 
derivative thus involves second derivatives of the velocity, requiring it to be 
at least of class C^. 

A second common type of derivative is the corotational or Jaumann derivative 
(See Refs. [9] and [10, p. 342], and references therein), where the local vortic- 
ity of the flow is incorporated into the derivative operator. The corotational 
derivative is compatible with the metric, but like the convective derivative it 
depends on gradients of the velocity. 

The third type of derivative we discuss is a new, time-dependent version of the 
usual directional derivative along a curve used to define parallel transport [3- 
5]. The curve here is the actual trajectory of a fluid particle, with tangent 
vector given by the Eulerian velocity field. The directional derivative does 
not depend on gradients of the velocity field. The concept of time-dependent 
parallel transport can be introduced using this derivative, and is equivalent to a 
covariant description of advection without stretching. A directional derivative 
was introduced in the context of fluid motion by Truesdell [11, p. 42], but it 
does not allow for time-dependence in the coordinates or metric. (Truesdell 
calls it the material derivative because of its connexion to fluid elements.) 

In this paper, we present a unified derivation of these different types of covari- 
ant time derivatives. We do not restrict ourselves to Eulerian and Lagrangian 
coordinates, as this obscures the general covariance of the theory: both these 
descriptions lack certain terms that vanish because of the special nature of 
the coordinates. From a dynamical system defined in some Eulerian frame, we 
transform to general time-dependent coordinates. We then find a transforma- 
tion law between two time-dependent frames with no explicit reference to the 
Eulerian coordinates. The Eulerian velocity of the flow is not a tensor, but 
the move to general coordinates allows the identification of a velocity tensor 
that transforms appropriately (Section 2). We also derive a time evolution 
equation for the Jacobian matrix of a coordinate transformation between two 
arbitrary time-independent frames. This time evolution equation facilitates 
the construction of the covariant time derivative in Section 3. After a discus- 
sion of the rate-of-strain tensor in Section 4, we present in Section 5 the three 
types of covariant time derivatives mentioned above: convective, corotational, 
and directional. 

Section 6 addresses a fundamental issue when dealing with generalised coor- 
dinates: the problem of commuting derivatives. In manipulating fluid equa- 
tions it is often necessary to commute the order of time and space derivation. 
When commuting two covariant spatial derivatives, the Riemann curvature 
tensor must be taken into account. Similarly, when commuting a covariant 
time derivative with a spatial derivative, there arises a tensor we call the 
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time- curvature. This tensor vanishes for sufficiently simple time-dependence 
of the metric, and satisfies many properties similar to the Riemann tensor. 

Throughout this paper, we will usually refer to the "fiuid," "fiuid elements," 
and "velocity," but this is merely a useful concretion. The methods developed 
apply to general dynamical systems where the velocity is some arbitrary vector 
field defined on a manifold. The covariant time derivative still refers to the rate 
of change of tensors along the trajectory, but the tensors do not necessarily 
correspond to identifiable physical quantities. For example, the covariant time 
derivative is useful in fornmulating methods for finding Lyapunov exponents 
on manifolds with nontrivial metrics [1]. 



2 Time-dependent Coordinates 



We consider the dynamical system on an n-dimensional smooth manifold U, 

x^v(t,x), (1) 

where the overdot indicates a time derivative and f is a differentiable vector 
field. (For simplicity, we restrict ourselves to a given chart.) A solution x{t) 
defines a curve C in U with tangent v. We view the x as special coordinates, 
called the Eulerian coordinates, and denote vectors expressed in the Eulerian 
coordinate basis {d/dx'^} by the indices i,j, k. 

A time-dependent coordinate change z{t,x) satisfies 

(2) 

where the d/dt\,j. is taken at constant x. Here and throughout the rest of the 
paper, we assume the usual Einstein convention of summing over repeated 
indices. We denote vectors expressed in the general coordinate basis {d/dz"'} 
by the indices a,b,c, d. We use the shorthand notation that the index on a vec- 
tor X characterises the components of that vector in the corresponding basis: 
thus X° and X* are the components of X in the bases {d/dz"'} and {d/dx^}, 
respectively. The components X" and X* are also understood to be functions 
of z and x, respectively, in addition to depending explicitly on time. 



Defining v :— z, we can regard Eq. (2) as a transformation law for v, 

(3) 



dz^ J. dz"' 



dx^ dt 

This last term prevents v from transforming like a tensor. (We refer the reader 
to standard texts in differential geometry for a more detailed discussion of 
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tensors [3-5].) 



Now consider a second coordinate system z{t,x), also defined in terms of x. 
We can use Eq. (3) and the chain rule to define a transformation law between z 
and z, 



(4) 



Any explicit reference to the coordinates x has disappeared (except in d/dz\x). 
Equation (4) is a transformation law between any two coordinate systems 
defined in terms of x, and implies that v"' — {dz°- / dt) \x transforms hke a tensor. 
This suggests defining the tensor 



dz" 


dz" 


1 -a dz^ 








dt 


X ~ dz'^ 


V dt 



dz° 



dt 



(5) 



which we call the velocity tensor. The velocity tensor is the absolute velocity 
of the fluid v with the velocity of the coordinates subtracted. 

In addition to the coordinates x, characterised by dx^/dt\x = 0, we intro- 
duce another special set of coordinates, the Lagrangian coordinates a, de- 
flned by d = 0. We denote vectors expressed in the Lagrangian coordinate 
basis {d/da'^} by the indices p and q. Prom Eq. (3), we have, 



0. 



(6) 



The initial conditions for a are chosen such that Eulerian and Lagrangian 
coordinates coincide at t = 0: a(0, x) = x. 

Lagrangian and Eulerian coordinates have the advantage that the time evolu- 
tion of their Jacobian matrix is easily obtained. The Jacobian matrix dx^da'^ 
satisfies [6] 



d ( dx^ 



dv^ dx^ 



By using the identity 



dt \ dai 



d ( dx^ daP^ 



dx^ da'i 



(7) 



dt \ da'i dx^ 



dt 



m = 0, 



which follows from the chain rule, and using the Leibniz property and Eq. (7), 
we find 



It Vdx^ 



da'^ dv^ 



dx^ dx^ 
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The Leibniz property can be used again to find the time evolution of the Jaco- 
bian matrix of two arbitrary time-dependent transformations z{t,x) and z{t,x), 

d I dz"-\ dv"" dz^ dz'^ dv^ 

'dt\d¥)~~d^'d¥~~d¥'d¥' ^' 

All reference to Eulerian and Lagrangian coordinates has disappeared from 
Eq. (8); this equation is crucial when deriving the covariant time derivative of 
Section 3. 



3 The Covciriant Time Derivative 



The standard time derivative operator, which we have been denoting by an 
overdot, is defined for a vector field X as 



dt 



b 



where recall that z — v. The first term is the change in X due to any exphcit 
time-dependence it might have; the second term is the change in X due to 
its dependence on z. The time derivative is not covariant, because a time- 
dependent change of basis will modify the form of Eq. (9). 

We define the covariant time derivative T) by 

M»:=X" + q;%X\ (10) 

where the q;"^ are time-dependent quantities that arc chosen to make DX" 
covariant. In order that the operator D have the Leibniz property, and that 
it reduce to the ordinary derivative (9) when acting on scalars, we require 

•DYa = Ya-a\Y,, 

when acting on a 1-form Y. When T) acts on mixed tensors of higher rank, 
an a must be added for each superscript, and one must be subtracted for 
each subscript. We refer to the a as connexions, by analogy with the spatial 
derivative case. 

By enforcing covariance of D, we can derive a general expression for a"'^- 
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Since X is a tensor, we can write 

M" = X> — X" 

d fdz^\ dz^ • ^ dz^ „ 

dz"- 
dz'' 

because is by definition covariant. Hence, we require the a's to transform 
as 

^ dz^ dz'' ^ dz" d I dz'^\ 

The first term in (11) is the usual tensorial transformation law. The second 
term implies that a is not a tensor, and arises because of the time-dependence. 

Inserting the evolution Eq. (8) into (11), we can rewrite the transformation 
law for a as 

= dv" dz" dz'' ( „ dv"\ 

implying that a"b + {dv" /dz'') transforms hke a tensor. Hence, 

dv" 

«^ = -^ + ^^, (13) 

where "K is an arbitrary tensor. Equation (13) is the most general form of the 
connexions a. 

In Section 5, we consider three convenient choices of the tensor "K. But first 
in Section 4 we examine the action of the covariant derivative on the metric 
tensor. 



4 The Rate-of-strain Tensor 



Our development so far has not made use of a metric tensor. We now introduce 
such a tensor, specifically a Riemannian metric g : Tli x TU The 
components gab of the metric are functions of z and t, but the indices a and b 
run over the dimension n of Tit, and so do not include a time component. 
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It is informative to consider the derivative of the metric tensor, 



"Dgab = gab - a" a Qh 

dgab 



dt 



+ 



OL hQac 
■V'' + 9ac^ + 9bc 



dz" 



i^ab + ^ba), 



where we have used the metric to lower the indices on IK. We define the 
intrinsic rate- of- strain or rate-of- deformation tensor 7 [8,7] as 



7a& = 



dg, 



ab 



dt 



(14) 



Here we denote by Va the covariant derivative with respect to 2;", 



dz^ 



The Riemann-Christoffel connexions are defined as [4,5] 



pa — 1^ 
^ be — 29 



1 ad ( ^dbd , dOcd dgbc 



dz' 



+ 



Q^b Q~,d 



(15) 



whence the identity 



i/acJ- bd i/focJ- ad ~ Q^d 



holds. The covariant time derivative of the metric can thus be rewritten 

'^Dgab = -fab-^lb, (16) 
where IH^ := ^{^ab + ^ba) denotes the symmetric part of "K. 

The rate-of-strain tensor 7 describes the stretching of fluid elements. The time 
derivative of the metric in its definition (14) is necessary for covariance under 

time-dependent transformations; the term describes straining motion that is 
inherent to the space, as embodied by the metric. The trace of the rate-of- 
strain tensor is a scalar 



7% = ^"'^7ac = Ve^^+2 



dt 



logbl , 



where \g\ is the determinant of gab and we have used the identity 

ac 9gac 
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dt 



d_ 
di 



The rate-of-strain tensor can be decomposed as 

lab = lab+^l\9ab, 



(17) 



(18) 



(19) 
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where is traceless and represents a straining motion without change of 
volume, and ^i^f-gab/n is an isotropic expansion. We see from the trace (17) 
that for a time-dependent metric there can be an isotropic expansion even 
for an incompressible flow, if d\g\/dt\z 7^ 0. Note also that in Lagrangian 
coordinates (characterised by v'^ — Q), the rate-of-strain tensor reduces to 

SO that the deformation of the space is contained entirely in the metric tensor. 



5 Three Covariant Derivatives 



As mentioned in Section 2, the requirement of covariance only fixes the covari- 
ant time derivative up to an arbitrary tensor [Eq. (13)]. That tensor may be 
chosen to suit the problem at hand, but there are three particular choices that 
merit special attention. In Section 5.1 we treat the convective derivative, and 
in Section 5.2 we examine two types of compatible derivatives: corotational 
and directional. 



5.1 The Convective Derivative 



The choice !K = is equivalent to the convective derivative of Oldroyd [6,7]. 
The connexion, Eq. (13), reduces to the simple form 



The convective derivative Vc acting on a vector X" is thus 



dt 



+ 



(20) 



When acting on a contravariant vector X", as in Eq. (20), !Dc is sometimes 
called the upper convected derivative [9]; acting on a covariant vector 1^, 
DcYa = Ya+ {dv^ /dz"-) Ffe, is then called the lower convected derivative. 



In general, for an arbitrary tensor T, 



d7_ 
'dt 
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Table 1 

Comparison of the equation of motion for the components of an advected and 
stretched vector field B. The equations for the covariant and contravariant com- 
ponents of VcB differ because of the lack of compatibility with the metric. 

Type Contravariant components Covariant components 

Convective 2)c5« = Dc5„ = 2Sc7^„ 

Corotational 'DjB°- = B''-fc'' 'DjBa = Bcja'' 

Directional 'D,B" = B'^ (V, V" + /-,:,") 'D,B„ = B, (V V; + h,/) 



where is the Lie derivative of T with respect to v [4,5]. In Lagrangian 
coordinates, we have = 0, so the convective derivative reduces to 



The convective derivative is not compatible with the metric: from Eq. (16), 
the metric's derivative is 

I'^cQab = Tab ■ 

which does not vanish, unless the velocity field is strain-free. 

The convective derivative is ideally suited to problems of advection with 
stretching, where a tensor is carried and stretched by a velocity field. Ta- 
ble 1 summarises the form of the equation for advection with stretching of a 
vector field B (B is "frozen in" the flow [12]) for the three different types of 
derivatives introduced here. The equation for the contravariant component B^ 
is simply Dc-B" = 0, but the equation for the covariant component Ba = Qac B^ 
is DcBa = 2Bc'j'^a- These two equations differ because the operator Dc is not 
compatible with the metric. 



5.2 Compatible Derivatives 



Another way to fix J-C is to require that the operator D be compatible with 
the metric, that is, D(?a6 = 0. This allows us to raise and lower indices through 
the operator D, a property possessed by the covariant spatial derivative. Prom 
Eq. (16), the requirement Dgab — uniquely specifies the symmetric part 
of "Kab, so that "K^ — 7. Using Eqs. (13) and (14), we then find 

- i [gacVbV' - Qbc^aV'] + , (21) 

z 



Oiah 



^2 Qt 
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where := ^{"Kab — '^iba) is the antisymmetric part of "K. 
We define the antisymmetric vorticity tensor 



- 9bc 



(22) 



and the symmetric coordinate rate-of-strain tensor 



(dz^ 




(23) 



In Eulerian coordinates, we have Kij — ^{dgij/dt)\x- The compatible connex- 
ion (21) can be rewritten 



' dz'^ 
'dt 



+ l^ab — <^ab + ■ 



(24) 



Since CK^^ is antisymmetric, we can use it to cancel the vorticity, or we can set 
it to zero. The two choices are discussed separately in Sections 5.2.1 and 5.2.2. 

The decomposition of the velocity gradient tensor W into the rate-of-strain 
and vorticity tensors has the form 



9ac ^bV = [Tab - Kab] + iOab 



(25) 



in general time-dependent coordinates. When the coordinates have no time 
dependence, the tensor k vanishes, as does the derivatives dz^/dt\xi and we 
recover the usual decomposition of the velocity gradient tensor into the rate- 
of-strain and the vorticity. We can think of k as the contribution to the rate- 
of-strain tensor that is due to coordinate deformation and not to gradients of 
the velocity field. However, the term dg / dt\z \s a, "real" effect representing the 
deformation due to a time-dependent metric, and is thus also included in the 
definition of the intrinsic rate-of-strain tensor, 7, defined by Eq. (14). 

In Euclidean space, when the ratc-of-strain tensor 7 vanishes everywhere we 
are left with rigid-body rotation at a constant rate given by u) [7]. With an 
arbitrary metric and time-dependent coordinates the situation is not so sim- 
ple: the very concept of rigid-body rotation is not well-defined. Hence, even 
when 7 = 0, we cannot expect to be able to solve for v in closed form. 



5.2.1 The Corotational Derivative 

In this instance we choose the antisymmetric part !K^j to be zero. We call 
corotational the resulting covariant derivative, and denote it by Dj (the sub- 
script J stands for Jaumann). The appellation "corotational" really applies 
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to the Euclidean limit, gij = Sij, for which the compatible connexion Eq. (21) 
reduces to a.;,- = —ujij. It is then clear that the covariant derivative is designed 
to include the effects of local rotation of the flow, as embodied by the vortic- 
ity. (See Refs. [9] and [10, p. 342], and references therein.) The derivative (21) 
with = is thus a generahsation of the corotational derivative to include 
the effect of time-dependent non-Euchdean coordinates. 

In Table 1, we can see that, written using Dj, the equation for advection 
with stretching of a vector B"' has the rate-of-strain tensor on the right-hand 
side. The "rotational" effects are included in !Dj, hence the terms that remain 
include only the strain. 



5.2.2 The Directional Derivative 

Another convenient choice is to set "K^^ = cUab, thus cancelling the vorticity in 
Eq. (24). The resulting covariant time derivative then has the property that, 
in the absence of any explicit time-dependence, it reduces to the covariant 
derivative along the curve C [4,5], or directional derivative, where C is the 
trajectory of the dynamical system in the general coordinates z (Section 2). 
The derivative is called directional because it only depends on v, and not 
gradients of v. 

The form of the equation for advection with stretching of a vector written 
using is shown in Table 1. The W term on the right-hand side is the 
"stretching" term [12] (called vortex stretching when B is the vorticity vec- 
tor [13]). The K term represents coordinate stretching, and does not appear in 
Euclidean space with time-independent coordinates. 

Because the directional derivative depends only on v and not its gradients, it 
can be used to define time-dependent parallel transport of tensors. A vector X 
is said to be parallel transported along v if it satisfies X = 0, or equivalently 



dt 



+ V^X" = X' 



c 



dt 



(26) 



This can be readily generalised to tensors of higher rank. In Euclidean space, 
with time-independent coordinates, the right-hand side of Eq. (26) vanishes, 
leaving only advection of the components of X. Thus, parallel transport is 
closely related to advection without stretching; Equation (26) is the covariant 
formulation of the passive advection equation. 
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6 Time-curvature 



A hallmark of generalised coordinates is the possibility of having nonzero cur- 
vature. The curvature reflects the lack of commutativity of covariant deriva- 
tives, and is tied to parallel transport of vectors along curves [4,5]. An anal- 
ogous curvature arises when we try to commute D and V, respectively the 
covariant time and space derivatives: 



+ 9 



he 



Va {Oicd — led — ^ed) + RedaeV^ + \ Scda X'^ , (27) 



where the time-curvature tensor is defined by 

dgeb 

-rgbe'^e I - 

'dge 



'abc 



dt 



- V, 



dt 



+ Qae^c 



'dt 



+ R, 



■abee 



dt 



and the Riemann curvature tensor R obeys [5] 



(28) 



(29) 



The timc-curvaturc tensor satisfies Sabc = —Shaei ^'^^ the Riemann curvature 

tensor satisfies Rabed = -Rbaed, Rabed = Redab- 

Even for trivial (Euclidean) coordinates, we do not expect D and V to com- 
mute, because of the derivatives of v in the Va {^cd — led — ^cd) term of Eq. 
(27). Note that the coordinate rate-of-strain tensor defined by Eq. (23), 
does not appear in Eq. (27). 

The VX term in Eq. (27) vanishes for the convective derivative of Section 5.1, 
since then J{ = 0. For the directional derivative of Section 5.2.2, we have 
= led + '^ed-i SO the Commutation relation simplifies to 



be 



RedaeV^ + o s, 



2 '-^cda 



which does not involve second derivatives of v. For the corotational derivative 
of Section 5.2.1, with Ji^i = led, no terms drop out. 

The terms involving !K in the commutation relation (27) reflect properties 
of the velocity field v. In contrast, the tensors R and S embody intrinsic 
properties of the metric tensor g. The Riemann tensor R is nonzero when the 
space is curved. The time-curvature tensor S is new and has characteristics 
that are analogous to the Riemann tensor. It satisfies a cyclic permutation 
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identity, 




(30) 



which corresponds to the first Bianchi identity of the Riemann tensor. The 
time-curvature does not appear to satisfy an analogue of the second Bianchi 
identity. 

The property Sabc = —Sbac, together with the Bianchi identity (30), imply that 
S has n{n^ — l)/3 independent components, compared to the n'^{n'^ — 1)/12 
components of -R, where n is the dimension of the space. Thus one-dimensional 
manifolds have vanishing S and R. For 1 < n < 3, 5* has more independent 
components than R; for n — A, they both have 20. For n > A, R has more 
independent components than S. 

The time-curvature S vanishes for a time-independent metric and coordinates. 
It also vanishes for a metric of the form gij{t,x) = P{t) hij{x), where /i is a 
time-independent metric and x are the Eulerian coordinates. It follows from 
its tcnsorial nature that the time-curvature must then vanish in any time- 
dependent coordinates. In general, it is convenient to find S in Eulerian coor- 
dinates (where dx/dt\x = 0), 



and then transform Sijk to arbitrary time-dependent coordinates using the 
tensorial law. 



7 Discussion 

In this paper, we aimed to provide a systematic framework to handle com- 
plicated time-dependent metrics and coordinate systems on manifolds. The 
explicit form of the relevant tensors is often fairly involved, but the advantage 
is that they can be evaluated in time-independent Eulerian coordinates and 
then transformed to arbitrary coordinate systems using the usual tensorial 
transformation laws. 

The covariance of the time derivatives is made explicit by using arbitrary 
time-dependent coordinates. The results for the Eulerian coordinates are 
recovered by setting dx^/dt\x — 0, and those for the Lagrangian coordinates a* 
by setting v"^ — 0. 

The introduction of the time-curvature tensor allows us to treat the temporal 
dependence of the metric tensor in a manner analogous to its spatial depen- 




(31) 
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dence. For simple time-dependence, the time-curvature vanishes, such as for 
the case of a time-independent metric multiphed by a time-dependent scalar. 
As for the (spatial) Riemann curvature tensor, the components of the time- 
curvature can be computed for a given metric, and then inserted whenever a 
temporal and spatial derivative need to be commuted. 

We have only addressed the kinematics of fluid motion. The dynamical equa- 
tions relating the rate of change of quantities to the forces in play have not 
been discussed (see Refs. [9,6,8,7]), and depend on the specifics of the problem 
at hand. Nevertheless, covariant time derivatives provide a powerful framework 
in which to formulate such dynamical equations. 
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